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We study the propagation of the light mesons <j,ui,p, and ao(980) in dense hadronic 
matter in an extended derivative scalar coupling model. Within the scheme proposed it is 
possible to unambiguously define effective density-dependent couplings at the lagrangian 
level. We first apply the model to study asymmetric nuclear matter with fixed isospin 
asymmetry, and then we pay particular attention to hypermatter in /3-equilibrium. The 
equation of state and the potential contribution to the symmetry coefficient arising from 
the mean field approximation are investigated. 
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■ I. INTRODUCTION 
PLh " 

, In recent years the medium dependence of the meson-baryon couplings has been object of speculation 

This subject has been promoted by the successes of the so-called quantum hadrodynamics theory 
(QHD) jlo|, in the relativistic description of diverse nuclear phenomena. 
^ ' The assumption of variable couplings in the mean field approximation (MFA) is founded on different 

\ grounds. It can be interpreted as the trace of the quark structure of hadrons iQ^, or it can be viewed 

■ as a way to match effective lagrangians and free-space nucleon-nucleon interactions In any case 
' assigning a variable behavior to the couplings seems to be an appropriate method to interpolate from one 

CSj , dynamical regime to another, using effective hadron field models. A similar meaning has been given to the 
in- medium meson masses, which have been related to the transition to the chiral regime The so-called 
, Brown-Rho scaling law qualitatively describes the behaviour of the hadronic masses in the proximity of 
' the transition point. According to the scaling law, all hadronic masses decrease approximately at the 
\ same rate as the system approaches to the chiral phase transition (with exception of the pseudo-scalar 
meson masses). Applied to the light vector mesons this hypothesis could explain some experimental 
O ' results, as for example the dilepton production rate in heavy ion collisions. 

^ . On the other hand, in certain purely hadronic models [ p^|JT^ ] the non-polynomial meson-nucleon inter- 

C ■ action gives rise to effective, density-dependent coupling in the MFA [^-16|. In this paper we propose an 
extension of the derivative scalar coupling model (DSCM) of which preserves the charge symmetry 
and provides effective couplings for all the mesonic channels. The medium dependence comes through the 
mean value of the scalar a meson, evaluated at finite density and temperature. More general interactions 
could include non-linear vertices in terms of all the mesons considered, but we restrict here to the sim- 
plest case. This choice is similar to the approach of where the functional dependence of the couplings 
includes only the product ^'7^^', but not the scalar 'J'^' or more involved nucleon field combinations. As 
stressed in this reference, the density variation of the couplings must be written as a Lorentz invariant 
functional of the considered fields, to obtain the correct Euler-Lagrange equations. Otherwise the so 
called "rearrangement" contribution is absent and the thermodynamical consistency is lost. 

We use the model proposed here to investigate asymmetric nuclear matter at finite temperature, there- 
fore we explicitly include the isovector mesons p and ao(980). 

The equation of state of asymmetric matter is an important input in astrophysical studies such as the 
cooling rate of neutron stars or the supernova collapse mechanism. Also it is of interest in the description 
of heavy-ion collisions, where experiences with radioactive beams are expected to provide new insights 
into the structure of dense matter with high degree of asymmetry. This situation has stimulated the 
interest in this issue, and many theoretical works have been done in the last years |^ 17-25]]. 

As an important case of asymmetric matter we consider hadronic matter in /^-equilibrium. For this 
purpose we generalize the DSCM model to include the octet of baryons n, A, S, and S. It is well known 
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that the proton concentration is determinant in the cooling of neutron stars 1 26 2j] , and this concentration 
is mainly determined by the isospin-dependent contribution to the equation of state. 

It is also of interest to study the meson properties under these conditions since, for instance, it could 
carry information about transient states of matter during heavy-ion collisions. This subject has been 
poorly developed in the literature. We have evaluated the meson propagators in the relativistic random 
phase approximation (RRPA), including particle-antiparticle contributions, and we have extracted from 
them the effective meson masses. 
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We organize this paper presenting the model in Sec. |l|, in Sec. 11] we discuss the bulk properties of 
symmetric nuclear matter at zero temperature, meanwhile the asymmetric nuclear matter equation of 
state is treated in Sec. [V. The properties of /3-stable matter are considered in Sec. ^ and the Feynman 
graphs contributing to the RRPA, the evaluation of the propagators and the behavior of mesons in the 
hadronic environment are presented in Sec. 



VI. We conclude with the discussion and summary in Sec. 



VII. 



II. THE MODIFIED DSCM 



In this section we present a relativistic model of hadronic fields inspired on the DSCM proposed by 
Zimanyi and Moszkowski The DSCM has been used to study nuclear many-body effects in several 
applications |p8||, to investigate neutron star properties |p9| , extended to include nucleon resonances 
and hyperons ]3l| , related to an effective quark description of hadronic properties | |3^ , and generalized 
with a tensor coupling in order to improve the spin-orbit splitting. 

The DSCM has two important features which distinguish it from the QHD-I model of Ref. In first 
place it is non-renormalizable ab initio and there is no immediate way to introduce vacuum corrections to 
the MFA to the ground state, although the main properties of nuclear matter are successfully described. 
In second place a residual interaction can be extracted beyond the lowest order solution, whose strength 
decreases monotonically as a function of the baryonic density |jl^,|l6| . This fact ensures the ground state 
predominance at high density as assumed in QHD ]Tot . 

Since we are interested here in the description of asymmetric matter, besides the fields Vt" for the 
nucleons, we include the isoscalar scalar (a) and isoscalar vector (w^) mesonic fields, and those corre- 
sponding to the p isovector vector (p^) and the ao(980) isovector scalar {S^) mesons. We use greek, latin 
lowercase and latin uppercase indexes to denote Lorentz, baryon isospin and meson isospin components, 
respectively. 



In its simplest version the DSCM 1 13 1 has a Yukawa type N-w coupling and a N-cr non-polynomic term. 
We modify the vertices allowing for two different mesons (one of them the scalar a) to locally interact 
with a baryon: 
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where '^{x) is the isospin multiplet nucleon field, M is the averaged nucleon mass and gs,gd,gv, and gr 
are adimensional coupling constants. As usual in QHD the ground state for homogeneous infinite matter 
is approximated by considering mesonic fields as classical quantities and assimilating them to effective 
nucleon properties. Thus we can separate the c-number contributions: 



(t{x) ~ a + s{x), 
S^{x)=Sd'^'^ + d^{x), 

UJf_,{x) = CJ (5^0 + w^{x), 

p^{x)^p5,oS'^ + r^{x), 
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and are quantum fluctuations which 



where ct, S, cj, and p are classical mean field values and s, , 

are not included in the ground state. Expressions for the c-number contribution to meson fields can be 
obtained by taking statistical averaged Euler-Lagrange equations, and requiring self-consistency. In this 
way we obtain: 
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where we have used = 1 -|- g^a/M . The expectation values must be evaluated with the ground state 
solution for the nucleon field, which depends on a and 6 through the effective nucleon mass: 



(10) 



with li — 1,-1 for protons and neutrons, respectively. The nuclcon dispersion relation is also modified 
according to (po — 9w^jj — grlip/'^Y ~ = M*"^. 

In Eqs. (|^) and the terms between angular brackets represent the conserved baryon density and 
the isospin density, respectively. 

A residual nucleon-meson interaction arises beyond the lowest order approximation [ p^ by inserting 
Eqs. (||)-(||) in the interaction term: 

l + gsCj/M iV(l + fs) 

In the expression above the symbol 4>i represents any one of the mesonic fields 5, and p, which by 
virtue of Eqs.(|^)-(|^) splits into the classical mean value 0i, and the quantum fluctuation 5(j}i. Ti stands 
for the bare meson nucleon vertices: Ti — —g^r, g^^, corresponding to the aQ,LU, and p mesons 

respectively. The right hand side of this equation is non-polynomic and can not be used to directly 
apply a diagrammatic expansion. Restricting to the physical regime for which quantum fluctuations are 
negligible compared to mean values, enables us to approximate: 

^i±I^o.M*+joSe + Cres, (11) 
1 + gscr/M 

with 

fc = ^r*0„ (12) 
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We have introduced the medium dependent vertices F*, which are obtained from the bare ones by 
replacing the coupling constants gd,gw, and gr by effective couplings. The last ones are given by the 
relation gd/ g*d — gw/ gtu — gr/ g*- — N. Also we have used gs/gt = N'^ 



The expansion proposed in Eq. (11) respects the organizational principle of nuclear effective field 
theories Q. 

In this approximation, the residual interaction of Eq. ( |l^ ) arises besides the nucleon effective mass (|o|) 
and the contribution to the nucleon single particle energy (|l2|). The interaction term Cres comprises a 
one meson-nucleon vertex, together with a two-meson exchange term. In all cases the vertex functions 
are medium dependent, 5* (l + J2i ^i^fi/M) , F^, F^, and F* for the one- meson case and g*Ti/M for the 
two-meson instance. 

This linearized version can be used to study the quantum corrections beyond the mean field approxi- 
mation. 

Variable couplings are an expected feature of hadronic models, whenever the quark substructure be- 
comes relevant ftfl. Furthermore, density dependent couplings have been proposed as a key assumption 
in order to match relativistic nucleon potentials adjusted to scattering data, with hadronic field models 
[^]-^. This approach was initiated as a way to avoid involved Brueckner-Hartree-Fock calculations for 
finite systems, using one boson exchange potentials. Thus, the main purpose is to take advantage of 
the relative simplicity of the Hartree approach to the QHD models. The link between both schemes is 
established by requiring the equality of the nucleon self-energy in symmetric nuclear matter as evaluated 
in both formulations, and allowing QHD coupling constants to be density dependent. 

In our treatment the effective couplings are unambiguously extracted from the lagrangian, once the 
MFA has been introduced. Thus in this scheme one has a well defined and invariant way to describe the 
medium infiuence on the couplings. Furthermore, the internal consistency of the approach is guaranteed. 

Up to this point we have restricted the discussion only to protons and neutrons, however the intro- 
duction of hyperons is straightforward. A sum over different baryonic species must be considered in the 
lagrangian density and the vertices must be modified in order to take into account the isospin degeneracy 
of each one. Also, additional couplings between the mesons and every hyperon iso-multiplet must be 
introduced. A more detailed discussion will be given in Sec. 0. 

The energy density E for infinite homogeneous hadronic matter, can be evaluated in the MFA by taking 
the statistical average of the energy momentum tensor: E =< T^^ >. The thermodynamical pressure 
P under the same conditions is obtained by averaging the trace of the spatial-spatial component of this 
tensor: P =^<TtT^^ > /3. We include the corresponding equations for the sake of completeness: 



q 



i—p,n 



d^kEk^[nF{Ek^) + nFi-Ek^)] + -{mla^ + mlS^-mlL^^-mlf) 
(^■^y Jo ^ 



+ 9w^ n + 25rP ("p - "n), 



. 1 /-oo .2 



(14) 
(15) 



with Ek^ = ^JM*^ + k^, M* is given by Eq. (|T^), and 
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is the nucleon statistical occupation number, Xi ~ Zi + gwu) + Iigrp/2 is the particle energy, /i,; is the 
chemical potential, and /3 = l/ksT. The chemical potential is related to the number density of the i-type 
particle through: 
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and finally n — Up + Un is the total particle number density. 

To evaluate Eqs. (H) and (p^, one must fix the particle number densities np,nn and then simultane- 
ously solve Eqs. (|6|)-(|1C1|), together with Eq. (p^). 

Another interesting quantity is the nuclear symmetry energy Eg defined as: 
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where x = (^n ~ ^p) /it-- This energy contains a purely kinetic term T, and a contribution Vg coming from 
the isovector mesons only. The explicit expression for Eg at zero temperature, including the uq meson 
contribution has been already derived, see for example pi[: 
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where is the Fermi momentum for the i-type particle, Ep^ = \ Pp. + M* , and 
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III. BULK PROPERTIES OF SYMMETRIC NUCLEAR MATTER AT T=0 



In the previous section we have presented the model, which contains several free parameters. The 
masses of the ao, w and p mesons are taken at their physical values = 984 MeV, = 783 MeV, and 
rrir — 770 MeV respectively. We adopt the accepted value for the a meson mass nig — 550 MeV. There 
remains to determine the four coupling constants. We adjust them to reproduce the main bulk properties 
of symmetric nuclear matter: the saturation density uq — 0.15 fm~^ , the binding energy eb = —15 
MeV, and the symmetry energy Eg = 32 MeV at zero temperature and at normal density. Another 
quantity of physical interest is the isothermal compressibility kt, however the DSCM provides very good 
values for kt without imposing any further condition. Therefore we have three physical conditions to 
fix the four coupling constants. Two of them, gg and g^ are univocally determined to take the values 
gg = 12.379, gw = 14.624, whereas gr and gd are functionally dependent through Eqs. (17)-(^9|) evaluated 
at n = ng. For our calculations we have selected two sets of couplings (gnSd); denoted by A and B: set 
A = (11.583, 0) and set B = (15, 6.538), which are shown in Fig. 
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FIG. 1. The relationship between the couplings Qr and ga, constrained to reproduce the symmetry energy 
Es = 32 MeV at zero temperature and at the saturation density no. The pair of couplings A and B used in our 
calculations are marked with circles. 

As previously mentioned, a feature of the model proposed is the presence of effective couplings. The 
behaviour of these couplings relative to their vacuum values is the same for the (5, and p fields, and 
different for the scalar cr, as discussed in Sec. ||. As can be seen in Fig. ^ the channel corresponding 
to the last case is much more suppressed in dense matter. This figure corresponds to symmetric nuclear 
matter, but it must be noted that the behavior of the effective couplings depends on the composition of 
the hadronic medium, i.e. they must depend on the asymmetry coefficient x- 

Although of diverse inspiration and derivation, we compare these results with the density dependent 
hadron field theory (DDHFT) outcomes For this purpose we use the interpolating algebraic function 

given in Differences are appreciable at medium and high densities. The couplings for cr,uj, and p 
are monotonous decreasing in both formalisms. A dropping of 20% for the isoscalar mesons, and of 42% 
for the p meson is detected in the DDHFT at n/no ^ 2. Our results provides for the same conditions a 
stronger decay of 60% for the a coupling and of 40% in the remaining cases. 
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FIG. 2. The effective couplings relative to their vacuum values in symmetric nuclear matter at T=0, in terms 
of the baryon density. The solid line corresponds to the a channel, and the dashed line to the S,u), and p cases. 
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FIG. 3. The binding energy (left) and the pressure (right) as functions of the baryon number density at T=0, 
in the MFA. Solid and dashed lines correspond to our results and the standard DSCM calculations, respectively. 

With the sets of parameters A and B we have evaluated some properties of symmetric nuclear matter at 
zero temperature. In Fig. ^ we compare our results for the binding energy and the pressure as functions 
of the baryon number density, with the corresponding outcome of the standard DSCM. It can be seen 
that there are not appreciable differences below n = 1.5 no, from here on both eb and P grow more 
slowly in our calculations. The isothermal compressibility is a measure of the stiffness of the pressure, we 
get at the saturation density kt — 165 MeV, against kt — 220 MeV for the DSCM. The lower slope of 
the binding energy in our results is essentially due to the weakening of the repulsion at higher densities 
induced by the normalization factor N . On the other hand, the relative difference between the mean 
values a and u) increases with n in our model, meanwhile in the standard DSCM it approaches to zero. 
This gives rise to the relative lessening of the pressure at high densities in our results. 

The medium effects on the effective nucleon mass can be seen in Fig. ^, where a comparison with the 
DSCM result is made. In both cases M* is positive definite and monotonous decreasing, but the rate 
of falling at densities < n < 2 no is more pronounced in our case because of the higher value of gs 
needed to reproduce the normal properties of nuclear matter. At higher values of the baryonic density 
M* stabilizes, due to the dynamical screening of the effective coupling. 
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FIG. 4. The effective nucleon mass at zero temperature in symmetric nuclear matter. The line convention is 
the same as in Fig. H 
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FIG. 5. The full symmetry energy (left panel) and the interaction contribution Vs (right panel), in terms of 
the baryon number density. In both cases solid lines correspond to results with coupling sets A and B. Dashed 
hues are used in the left panel for the QHD-I model (without 5 contribution), and the parameterizations PI, 
P2 and P3 as explained in the text. In the right panel dashed lines correspond to the QHD-I model and to the 
parameterization given in Ref. for variational calculations (VC). 

As the next step we investigate the density dependence of the symmetry energy Eg. It has been 
profusely studied in the past, using non-relativistic potentials as well as relativistic formulations |p"7|-^3||. 
Recently Eg has received attention by its applications to the study of the structure of nuclei with a 
large neutron excess, produced in the radioactive ion beam facilities. It is also a relevant subject in the 
evolution of neutron stars, determining the composition of the ground state and the cooling mechanism 
[p6| , or the phase transition to quark matter [ p2| . Furthermore, it has been proposed that the ratio of 
neutrons to protons in the prc-equilibrium stage of collisions between neutron rich nuclei could distinguish 
the asymmetric contribution of the nuclear equation of state [ pT[ . 

Different theoretical predictions for Eg produce rather dissimilar density dependences. In the left part 
of Fig. ^we compare our results for the symmetry energy coefficient, with and without the contribution of 
the OQ-meson, with other commonly used descriptions. We include the result from the QHD-I model [0 
without scalar-isovector coupling, and other three cases labeled PI, P2, and P3. The latter correspond 
to the phenomenological parameterizations ||27||: 
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where u = n/no, Ci? is the non-relativistic Fermi energy at the saturation density, and the function F(u) 
takes the forms Fx(u) — 21*^/(1 -f u), F^iu) — u, and F3(u) — y/u for the curves denoted as PI, P2 
and P3, respectively. All the curves are almost coincident for densities n < 1.5 np, but their mutual 
differences become significant for densities above that limit. The exception corresponds to the cases A 
and P3, which differ each other only by negligible amounts in all the range of shown densities. From Eq. 
( p!9| ) it can be seen that the contributions to Vg of the iso-vector S and p mesons are opposite in sign. 
However choosing gd brings on an enhanced behaviour of Eg, because the value of gr required to 
adjust Eg = 32 MeV at n = ng is bigger than gd (see Fig. p. It must be noted that the rate of growth of 
the cases A, B, and P3 decreases with density, whereas it remains approximately constant for the curves 
PI, P2 and QHD-I. 

The effect of polynomial self-interactions of the a field in QHD models has been studied in |2^. Both, 
the inclusion of exchange terms and of the 6 coupling in Hartree approximation enhance the density 
dependence of Eg , and therefore diverge from curves A and B in the left part of Fig. ^, being more alike 
to the PI parametrization. 

The behavior of Eg depends on the method of evaluation and the model of interaction used, the latter 
defines the Vg term. In the right part of Fig. |^ we display the interaction contributions for Eg obtained 
with the sets A and B, and we compare them with Vg extracted from the QHD-1 model and with 
the parameterization given by for the variational calculations (VC) made in It can be seen that 
our results are intermediate between QHD-1 and the VC results. A characteristic behaviour of the VC 



-7 



is that Vs becomes negative for densities bigger than certain typical value, causing the disappearance of 
protons in neutron stars at high densities. 

From the behavior of the symmetry terms shown in Fig. ^ we expect that the fraction of protons in 
star matter should be lower in our results as compared, for instance, with the QHD-I model prediction, 
although this fraction remains non vanishing for all densities in our case. The inclusion of the S coupling 
(curve B) slightly increases the presence of protons. 

IV. THE EQUATION OF STATE OF ASYMMETRIC NUCLEAR MATTER 

We study here the properties of nuclear matter at finite temperature by taking the asymmetry coefficient 
X as a free parameter. In the next section the isospin asymmetry will be determined by the conditions 
of electric charge neutrality and matter stability against electroweak decay. 

In first place we inspect the density dependence of the nucleon effective mass for fixed x- In Fig- |§ 
we compare results with and without S coupling at T=0 and x — 0.5. For gd = (set A) proton and 
neutron masses are degenerate, and for gd 7^ (set B) the neutron (proton) mass is lowered (enhanced) 
due to medium effects. The splitting is heightened as the density increases. 

Temperature effects are minimal in the range < T < 100 MeV, and more noticeable at high densities. 
For example, when the coupling set A is chosen, an increment of about 5 MeV in the nucleon effective 
mass is observed at n = 5no as the temperature is raised from T = to T = 100 MeV, at a given x- Of 
the same magnitude but opposite in sign is the effect of increasing the asymmetry from x = to x = 1 
at a fixed temperature. When the coupling set B is used, it is found that the in-medium mass splitting 
AM* = M* — M* decreases when the temperature is raised at fixed x- On the other hand, AM* is 
enhanced when the asymmetry is isothermally increased. Numerical values of this mass splitting depend 
on the set of couplings used, and we estimate the magnitude of both temperature and asymmetry effects 
calculating AM* with the set B at n = 5no. In this case the splitting reduces about 5 MeV in neutron 
matter when the temperature covers the range < T < 100 MeV, but an increment of approximately 50 
MeV is found in AM* if x is varied between x = snd x = 1 at fixed temperature. 

The thermodynamical pressure P has been evaluated using Eq. (p^), for several temperatures < 
T < 100 MeV, and several asymmetries < x < 1. As expected, increasing the temperature produces 
an enhancement of the pressure. This effect is strengthened by raising the asymmetry. The quantitative 
behavior of the pressure can be seen in Figs. Q and ||. In the first one we plot the pressure as a function 
of the number density at fixed asymmetry x = 0.25 and for several temperatures. For T > 20 MeV it is 
a monotonous increasing function of the density, whereas for T = it exhibits a region of instability for 
densities below uq. This instability gives rise to a liquid-gas phase transition |jl^. The results in Fig. |^ 
correspond to the set A. By using the set B qualitatively similar results are obtained. 
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FIG. 6. The efltective masses of proton and neutron as functions of the relative baryon number density for 
asymmetric nuclear matter at zero temperature. The solid (dashed) line corresponds to results with the coupling 
set A (B). 
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FIG. 7. The pressure P in terms of the nucleon number density at x = 0.25. The different curves correspond 
to temperatures T = 0, 20, 40, 60, 80, and 100 MeV. The arrow indicates the sense of growing temperatures. The 
results shown are obtained with the set A. 

The relevance of the asymmetry in our calculations can be observed in Fig. ^j. The higher the values 
of the asymmetry the stiffer the pressure raises, this effect being emphasized when the coupling gd is non 
zero. The liquid-gas instability remains for low T and n, disappearing for x close to 1. 
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FIG. 8. The pressure as function of the nucleon number density. Different types of lines correspond to the 
asymmetry coefficients X = 0, 0.25, 0.50, 0.75, and 1.0, respectively. The arrows indicate the direction of increasing 
X- The top (bottom) panel represents calculations for T=0 (T=100 MeV), and the left (right) side is devoted to 
results using the coupling set A (B). 
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V. HADRONIC MATTER IN /3-EQUILIBRIUM. 



The conditions in the interior of certain stellar objects like protoneutron stars, requires additional 
degrees of freedom to be included in the lagrangian density of Sec. ||. Due to the large densities reached in 
such systems, several physical phenomena could take place. The appearance of mesons and baryons with 
strangeness, pion and/or kaon condensation, the chiral symmetry restoration, and the phase transition to 
a quark-gluon plasma are some of the expected processes. They must be taken into account, in order to 
properly describe the high density behavior of the equation of state. In this section we complete the model 
proposed by including the hyperons A, S, and S, but we do not treat explicitly the chiral symmetry and 
quark degrees of freedom. Therefore our results should be valid until fluctuations preceding any phase 
transition become relevant. However we present here calculations in the range < n/riQ < 10 for the 
sake of comparison. 

As anticipated in Sec. ||, the modifications in the lagrangian density are straightforward since we 
retain the form of the interaction for all the baryons. A sum over the full octet N,A, E, and 2 must be 
considered in Eq. (|l|), and new couplings gs, gd, gw, and gr are introduced for the hyperons. Furthermore, 
the vertex between the p meson and the baryon B must be modified by including an appropriate coefficient: 
Ib3 = 1/2 for proton and S", Ibs — —1/2 for neutron and S~, Ib3 = 1 for S+, Ib3 = for A and S", 
and Ib3 = — 1 for S^. 

The new couplings should be fixed to reproduce some relevant quantity, according to the phenomeno- 
logical approach. We proceed in this way to determine the a- and lo-A couplings. Using hypernuclei 
data the A binding energy can be extrapolated to be £a = — 28 MeV at uq, thus we obtain gsA = 2.335, 
gwA = 2.099. For the other hyperons there are not accurate experimental data. Different arguments are 
commonly used to get numerical values, like SU{6) symmetry or vector meson dominance. For simplicity 
and to carry out computations, we adopt gs.wS = gs,wE = 5s,iuA and gr.ds = gr,d3 = 5r,dA = 5r,d, with- 
out any further justification. With this choice we obtain at uq very similar binding energies: es = —28.11 
MeV, and Se -28.27 MeV for the E and S hyperons. 

Neutron star matter is electrically neutral, by additional contributions coming from electrons and 
muons. Leptons are included by means of Dirac free particle terms in the lagrangian of Eq. (Q). The 
equilibrium for (3 decay imposes constraints among the baryon and lepton chemical potentials: fiB — 
IJ-n — 9s Me- Here we have used qb for the baryon electric charge in units of the positron charge, /ig, 
and p,B represents the chemical potentials for neutron, electron and the baryon B, respectively. On the 
other hand, the electric charge neutrality imposes = — J^i '^i + Ib^b, with ni and ub indicating 
the number density for leptons and baryons. At zero temperature we consider the Fermi momentum pp^ 
writing: m = p^pjiiir"^), m = y^Ppi + m'f for leptons, and = \/ppB + ^^b _+ 9*wB^ + 9rlB3P for 
baryons. The efective mass is a generalization of eq. (p^), Mg = [Mb — lB9dS)/NB with Ib ~ I for 
p, S+, and 5°, = for A and S", and Ib = —1 for n, E^, and S^. We have used Nb = 1 + gsB<^/MB, 
where Mb is the averaged mass of the baryon isomultiplet B. 




FIG. 9. The effective baryon masses in terms of the total baryon number density. The left (right) panel 
corresponds to calculations with set A (B). Different type of lines are assigned to each iso-multiplet, as indicated 
in both panels. 
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The effective baryon masses as functions of the baryonic number density are shown in Fig. In 
the results corresponding to set A, each isomultiplet remains degenerate in mass. The variation of the 
hyperon masses are much more moderate than for the nucleon masses. As a consequence of the specific 
interaction used, the heavier the baryon considered, the weaker the density dependence of its effective 
mass is. Using the couphng set B the isospin degeneracy is removed, enhancing or dropping the mass of 
particles with positive or negative isospin projection, respectively, as it is depicted on the right part of 

In Fig. |10|the relative population of the baryonic species is shown in terms of the total particle number, 
at T = 0. In the range of densities studied the full baryon octet is present, with exception of the S° 
when the coupling set B is used. The results obtained with gd = and gd = 6.538 are very similar for 
the leptons and the lightest baryons (p, n, and A). Differences between them become noticeable in the 
growth of populations of the heavier fermions E and S. The more obvious is the early appearance, at 
n ~ 3.3 no, and predominance of Y,~ particles in the results with the coupling set A. 

The equilibrium baryonic population is not perturbed by the presence of the 5 meson (set B) at low 
and medium densities. The effect of turning on the 5 interaction, is twofold and is emphasized at high 
densities. In first place the baryon-(5 interaction enhances the effective mass of S° and diminishes that 
of E~ and S~, increasing and lowering the corresponding thresholds. The more evident consequence of 
this is the absence of "EP particles in the range < n/riQ < 10 (right panel of Fig. [To| ). In second place 
the coupling g^ grows with gd, affecting more strongly to the iso-triplet S than the iso-duplet H, due to 
the factor 7^3. Since the p meson contribution to the chemical potential is positive and greater for 
than for S~, this causes the appearance of E^ and of S~ to be delayed and anticipated, respectively, 
going from the left to the right panel of Fig. The A and E° baryons, which do not couple to the 
isovector mesons, do not show appreciable changes in their distributions. On the other hand, from the 
comparatively earlier raising of the E"*" and S~ population obtained with set B, it is possible to infer 
that, in absolute values, the 5 contribution to the baryonic chemical potentials lies between one half and 
the total p contribution. Of course, these results are partially a consequence of our assumption of equal 
couplings gd and gr for all the hyperons considered. 




FIG. 10. Tlie fraction of leptons and baryon species present in hadronic matter in /3-equiIibrium at T = 0. 
The different curves are labeled with the corresponding particle name. The left (right) panel corresponds to 
calculations with the coupling set A (B). 

The pressure in terms of the baryon number density is exhibited in Fig. There are abrupt changes 
of slope in the curve corresponding to the set A, which coincide with the appearance of hyperons. Similar 
changes, but more attenuated, take place in the curve with the set B. 
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FIG. 11. The pressure of hadronic matter in /9-equilibrium at T = 0. As indicated in the figure, solid (dashed) 
line corresponds to calculations using the coupling set A (B). 



VI. MESON PROPAGATION IN ASYMMETRIC NUCLEAR MATTER 



Medium effects in the meson properties have received attention in the later years, as they could carry 
the signals of phase transitions in the hadronic environment. As previously stressed, we expect our results 
to be valid out of the vicinity of the transition point. 

In the MFA mesons are treated as classical fields, with constant mass. In order to include quantum 
corrections we must go beyond the MFA. This can be done in the relativistic random phase approximation 
(RRPA), using the linearized residual interaction of Eq. (p^). In this approach the meson propagators 
are corrected by incorporating the baryon bubble diagrams at all orders, by using the Dyson-Schwinger 
equation. From the corrected propagator the effective meson mass can be extracted. This procedure has 
been applied in QHD calculations, see fo r example ||l^ and references listed therein. Specific computations 
with the DSCM can be found in 



(a) 



(b) 





FIG. 12. Feynman diagrams included in the RRPA. Case a (b) corresponds to pure (mixing) meson propagation. 
The solid fine stands for baryon propagator, dashed and dotted lines represent meson propagators of different 
types, and the fiUed circle and the cross tfieir respective vertices. 



At second order, the one-loop proper polarization insertions comprise the diagrams shown in Fig. 
The case (a) represent the propagation of a pure meson field, and the case (b) the mixing amplitude of 
different mesonic types. Due to baryon current conservation, the proper polarization for all the mesons 
can be written in terms of a few components. Some of them are divergent and requires an appropriate 
regularization. For this purpose, we follow the scheme outlined in ||l^. We do not deduce those results 
here, but we reproduce the main equations for the sake of completeness. 

The formalism is best described within a generalized meson propagator, in a matrix representation of 
dimension equal to the sum of the mesonic degrees of freedom. For example the generalized free meson 
propagator V°, has in its diagonal blocks the free meson propagators S°{q), D^^^{q), W°y(q), and 
i?^^^(g), for the cr(x), (5(x), u;(x), and p{x) fields, respectively, and null matrices in the complementary 
spaces: 



/ S°{q) i 



KM 



D^{q) i 



W°{q) i 



1 o 



A similar expression holds for the full generalized propagator V , but the complementary spaces are 
filled with the mixing meson propagators: 



I S{q) M„s{q) M^^{q) M„p{q) 

Msa{q) D{q) MsUq) Msp{q) 

M^^iq) M^s{q) W{q) M^p{q) 

\Mp„{q) Mps{q) MpM) R{q) 



Q/3 



The Dyson-Schwinger equation can be used to solve for V ^(g): 

Kpi<l)^Kp\l)-'^c.p{q), 
where we have introduced the generalized polarization insertion 



na/3(g) = 



n,{q) n,s{q) TiaUq) 

Hs^iq) Ild{q) UsUq) ^Sp{q) 

^^a{q) ^w{q) M^piq) 

Hp^iq) Upsiq) UpUq) n,(g) 



a/3 



Since we are primarily interested in the propagation of the pure meson fields, we do not consider mixing 
polarizations. The formulae for the one-loop diagonal components are as follows: 



lU,{q) 



(2^ 



Tr [GBiq)GBiq+p)] 



J2^Tr[GBiq)T,GB{q 
^ Mb 



-P)] 



A, A' 



'A -PA 



Tr [GB{q)VxGB{q+p)Tx'] /Ml ) , 



mi^{q)^gfJ2 I [GB{q)T^GB{q + p)T^] , 

iIi':,%q)^Y.9ll j ^,Tr[GB{q)rGB{q+p)Y], 
il^^^^Aq) = 9;^Y1 1 704^^ [GB{qh,T^GB{q+p)rT^], 



where the index B runs over all the baryons considered, and A, A' in the first equation runs over the 
meson fields S,uj, and p. The vertices have been described in Sec. ||. The baryon propagators 
GB^q) are evaluated in the MFA. In our calculations we only need the transversal component in the 
Lorentz indices, and the third component of isospin. For this purpose we use = for the nucleon 
and E, — 1 for the A, and =diag(l, 0, — 1) for the S particle. The referred expressions contain 
particle-antiparticle, particle-hole, and Pauli blocking contributions. The former one is divergent, to 
extract finite contributions we apply the regularization scheme outlined in ||l5[] . The Lorentz scalar 
contributions, containing the integrand Tr [GBiq)GBiq + remains undefined by a constant A, related 
to the covariant derivative of the polarization evaluated at the regularization point. We take this constant 
as a free parameter to analyze the possible dynamical regimes. There are two independent parameters As 
and Ad corresponding to the diagonal components lis and 11^ respectively. We require null contribution 
for the polarization evaluated on the meson mass shell, at zero baryon density and temperature. Thus 
we obtain for the finite particle-antiparticle contribution of the baryon-B bubble: 



n' 33 



(q) 



27r2 

_2 



dz z{l 



In 



M^ 



Ml 



z(l-z)q2 



z(l — z)m^ 



KB{q) = ^c^(.rn*iml-e) 



35, 



*2 



47r2 



dz [M*^ 



z(l - z)q2] In 



z(l - z)q 



21 



Ml - z{l - z)ml 



(20) 
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FIG. 13. The effective a-meson mass as a function of tlie baryonic density at T=0. In tfie left (right) panel 
the results for the regularization parameter Xs = 10 (As = 100) are plotted. In each case the solid hne corre- 
sponds to hadronic matter in /3-equiIibrium, and the dashed lines to nuclear matter with asymmetry coefficients 
X = 0, 0.25, 0.50, 0.75, and 1. The arrow indicates the direction of growing x in the last case. 



where the index v = w,r runs over the vector mesons, and c = s,d runs over the scalar ones. In the case 

of isovector polarizations, it must be regarded as the (3,3) isospin component. Furthermore we have used 
= f/^q^, q is the modulus of the spatial component of the momentum, and to|j = Mg/Ms- 
Once the polarization has been properly defined, we introduce the effective meson masses m*, mjj, mj^, 

and m*. They have been defined as the zeroes of the corresponding inverse propagators at zero vector 

momentum, i. e. the pq solutions of: 

^aa'(PO,P - 0) = Vl~\po,P - 0) - n,,(po,P - 0) = 0, 

for a = s,d, w, and r. 

In Figs. [l3[|r^ we display the numerical results for the density dependence of the meson masses at 
T — 0, under different compositions of the hadronic medium. In Fig. |l3| the behavior of the a meson 
mass is presented. The results correspond to the coupling set A, there are no appreciable differences with 
respect to the calculations using set B . 




n/rio n/no 

FIG. 14. The effective ao-meson mass as a function of the density at T=0. The left (right) panel corresponds 
to the values — 10 (Ad — 100). Solid line stands for /3-stable hadronic matter, and dashed lines correspond 
to asymmetric nuclear matter with asymmetry x = 0.25,0.50,0.75, and 1. The arrows indicate the direction of 
growing x- 
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It can be seen that in /3-stable hadronic matter the mass is ahuost constant at high densities. In the case 
of nuclear matter at constant asymmetry the density dependence is more pronounced, and monotonous 
decreasing. The asymmetry dependence is smaU, as this figure shows. 

Fig. |l^ is devoted to the qq meson mass. The curves for matter in /3-equiHbrium are quahtatively 
similar to those corresponding to the a meson. In asymmetric nuclear matter its behavior is much more 
striking. For = 10 the mass becomes zero at n/no ~ 4, whereas for Xd = 100 it decreases sharply but 
never vanishes in the whole range considered. 

The masses for the vector mesons are shown in Fig. |l^. For the cu meson the behavior is almost 
independent of the composition of the hadronic environment, sensible departures are observed only for 
extreme densities. This is not the case of the p meson mass, a clear difference among /3-stable matter 
and asymmetric nuclear matter is shown, even at low densities. 

Since Ilaa{p) receives the contribution of all the baryonic species considered, the mesonic effective 
masses are strongly influenced by the inclusion of hyperons, even at densities close to the normal saturation 
value. It must be noted that the particle-antiparticle term coming from the baryon-B bubble contributes 
even when this particle is not present on its Fermi shell. The hyperon particle-antiparticle contributions 
at medium and high densities causes the stabilization of the mesonic masses in neutral /3-stable matter. 
This can be appreciated in Figs. |l3[]l5| where they exhibit a weaker density dependence as compared to 
nuclear matter results. The magnitude of this effect is distinct for each type of meson, depending on the 
strength of its coupling to the hyperons. This fact explains why the omega-meson mass is little affected 
by the presence of hyperons relative to the rho-meson mass. 
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FIG. 15. The density dependence of the vector mesons. The left panel corresponds to the iso-scalar lu meson, the 
right panel to the iso- vector p meson. In both cases the bold (thin) solid line corresponds to /3-stable matter with 
coupling set B (A), and dashed (dotted) lines correspond to asymmetric nuclear matter with x = 0, 0.25, 0.50, 0.75, 
and 1 using the set B (A). The arrow indicates the direction of growing x- 



In order to compare with other similar works, one can appreciate that at high densities m* is 
monotonous decreasing and remains below its vacuum value (see Fig. |l3|). On the contrary, the re- 
sults found in QHD-I model calculations for symmetric nuclear matter give increasing values of it for 
n > no 1^^. Meanwhile, the behavior of is more alike to ours. The inclusion of polynomic self- 
interactions of the a and uj fields [3^ ] produces increasing or fairly constant density dependence for , 
but if the direct a-uj interaction is omitted it results in a strongly decreasing behavior |37| . On the other 
hand, the p meson mass can be compared with the results shown in [Q. In both cases m* is below its 
vacuum value, for instance at the saturation density m* is about 10% down in our calculations against 
a decrease of 30 — 40% in [Q. The differences can be surely assigned to the absence of tensor couplings 
in our model. 

A comparison with the DSCM calculations of [|l5| shows that the main effect of generalizing the 
nonlinear a interaction to the vector meson couplings is noteworthy for the effective a and to meson 
masses. Their density dependence is more abrupt in the present work, and generally speaking the values 
obtained here are lower than in Ref lllSl. 
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VII. DISCUSSION AND SUMMARY 



In this paper we have proposed an effective relativistic hadronic model inspired in the DSCM to investi- 
gate in-medium hadronic properties, in terms of the baryon isospin asymmetry. The non-hnear cr-nucleon 
interaction is generahzed to the isoscalar vector, isovector scalar, and isovector vector channels. Effec- 
tive medium dependent couplings arise at the MFA, and a residual interaction with one and two meson 
exchange is obtained beyond the MFA. The equation of state (EoS) for symmetric nuclear matter is 
softer than the one corresponding to the DSCM. The symmetry energy coefficient shows an intermediate 
behavior between the QHD model and non-relativistic variational calculations. The asymmetry depen- 
dence of the EoS becomes relevant for densities n > 3no, and it is emphasized by the contribution of the 
ao(980) meson exchange. Temperature effects in the range < T < 100 MeV are noticeable in the EoS, 
but moderate in the effective baryon masses. As a particular manifestation of asymmetric matter we 
study hadronic matter with hyperons, in equilibrium against electroweak decay at T=0. The 6 coupling 
is the cause of notable modifications in the population of hyperons at high densities. Due to the lower 
hyperon-meson couplings, relative to the nucleon-meson ones, the hyperonic effective masses decrease 
more moderately as the baryon density increases. 

The effective meson masses have been evaluated at T=0 in the RRPA, including particle-antiparticle 
finite contributions. The regularization procedure left undefined parameters As and A^. We have selected 
numerical values for them, which differ by one order of magnitude, and are representative of the pos- 
sible dynamical regimes. In asymmetric nuclear matter the scalar a and 6 mesons exhibit monotonous 
decreasing masses for high densities, whereas the vector uj- and p-meson masses show a slight increase 
for n/riQ > 1. In all cases the effective masses remain below its vacuum values at extreme densities. The 
dependence on the asymmetry x is more evident for the vector mesons. In /3-stable hadronic matter the 
density variation of the effective masses of all the mesons considered is damped, becoming approximately 
constants at twice the saturation density. This effect is due to the particle-antiparticle terms, contributing 
even for particles out of their Fermi shell. 
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